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ABSTRACT 

Vertical plane response of submersible vehicles in the proximity of a free 
surface in both deep and shallow waters is evaluated using a potential flow, strip theory 
solver. Three criteria, namely periscope submergence, sail broaching, and collision are used 
to quantify the response. These criteria combined with the vehicle’s response amplitude 
operators in regular sinusoidal waves along with a statistical description of the seaway lead 
to an assessment of an overall operability index for the vehicle. The operability index is 
calculated within a given range for sea states and sea directions and for various vehicle 
speeds and operating depths. The results indicate that a certain combination of depth and 
speed can lead to a significant improvement in vehicle operations. 
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I. INTRODUCTION 



When a submarine is at periscope depth beneath a seaway, several exciting forces 
and moments are encountered. These include first-order oscillatory motions at the 
prevailing wave frequencies and second-order drifting motions at very low frequencies 
well outside the wave spectrum of the seaway. These are referred to as the free surface 
suction forces and moments. In practice these low frequency motions are difficult to 
control and may give rise to unsatisfactory depth keeping. Standard ways of computing 
free surface effects rely on combinations of potential flow and semi-empirical coefficient 
based models. Rankine type sources are distributed along the hull of the ship, which 
satisfy the free surface boundary condition. Source strength can be computed by 
satisfying the exact body boundary condition, that no fluid can pass through the hull 
surface. Discretization of the hull form into a finite set of Hess-Smith type quadrilateral 
panels allows formulation of algebraic system of equations to be solved for the unknown 
singularity strengths. [Ref. 1] Combination of forces and moments generated in this way 
with deep water force predictions can then be utilized to simulate the motion of the boat 
under waves. 

Although such suction forces and moments are difficult to predict, they are slowly 
varying in time and as a result they can be controlled by either the operators or automatic 
control systems. Techniques of disturbance estimation and compensation can be 
employed in order to allow satisfactory real-time estimation of unknown disturbances 
from depth, pitch angle, and pitch rate measurements [Ref. 2]. In contrast to these second- 
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order forces, first-order forces cannot be actively controlled in practice since they occur at 
very high frequencies (in the order of magnitude of a few seconds) which are normally 
outside the range of hardware response times. Therefore, it is essential that we have a 
clear understanding of the effects of first-order excitation on boat missions in order to 
maximize its window of operations. This approach will identify parameter regions where 
operations may be carried out with higher degrees of confidence of success and with 
better opportunities for active control. In order to achieve this first we need to have a 
computational tool which will allow us to determine vehicle motions in the vicinity of a 
free surface in deep and shallow water. In this work we utilize a strip theory seakeeping 
prediction program based on the work by Beck and Troesch [Ref. 3]. 

Accurate submarine maneuvering predictions are essential for operation, to 
provide optimal and safe submerged operating envelopes. Relative vertical motion is one 
of the most important response elements that appear to be the most crucial to the 
submarine operators for successful completion of a mission near a free surface. To decide 
whether the submarine could complete her mission successfully or not we need to adopt a 
number of criteria each pertaining to different operational hazards. These criteria can be 
broadly divided into two major categories, subtle and catastrophic failures. In this study 
we consider two criteria. One of them is periscope submergence, which is a subtle failure 
and the other is sail broaching, which is a catastrophic failure. 

Subtle failures refer to events which will occur in all types of periscope depth 
operations, i.e., propeller emergence, mast emergence, periscope submergence. Single 
occurrence of these events does not constitute failure of operations. However, their 



2 



frequency imposes operability limits for a certain sea state. Periscope submergence 
impairs visual information. The dominant criterion is the number of occurrences per unit 
time. 

Catastrophic failures will probably result in either cancel of operations, failure to 
complete mission or submarine detection, i.e., broaching, loss of depth control deep. 
Broaching the sail is defined here as any portion of the sail breaking the surface. It is 
assumed that submarine detection will occur with probability of one each time a 
broaching occurs. 

In shallow water operations, we employ an additional criterion in addition to the 
two above. This corresponds to the event of vehicle collision with the sea-bed and is 
determined statistically based on the absolute vertical motion of either the bow or the 
stem of the vehicle and its clearance from the sea-bed. Combination of all of the above 
criteria provides a quantifiable measure of vehicle operability within a certain range of 
sea-states and directions. 

Chapter II of this thesis describes the mathematical foundation for the problem. 
Chapter HI presents results for the operability index of a vehicle for periscope 
submergence and sail broaching criteria for different sea-states and wave heading angles 
and deep water operations. In Chapter IV the operability index of the vehicle is presented 
for periscope submergence, sail broaching, and collision criteria for shallow water 
operations. Finally, conclusions from this work and recommendations for further research 
are outlined in Chapter V. 
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II. EVALUATION OF RESPONSE 



A. GEOMETRY OF A SUBMARINE 

Froude (1877) [Ref. 4] introduced the concept of a ship with a forward end called 
the “entrance”, a parallel middle body, and an after end called the “run”. Chapman (1768) 
[Ref. 4] introduced the concept of a ship hull with the entrance a portion of an ellipsoid of 
revolution, and with the run a portion of a parabola of revolution. This concept is tailor- 
made for use in calculating volumes of modem submarine hulls, as described by Jackson 
(1983) [Ref. 4]. It was developed by assuming a body of revolution with a 
length/diameter {LID) ratio of six and a maximum diameter at 0.4L. The entrance has a 
length , L r , of 2.4 diameters. The run or after end has a length, L „ , of 3.6 diameters. The 
entrance can be calculated as an ellipsoid of revolution, and the run as a paraboloid of 
revolution which is rotated about a line parallel to the center-line. The equations of the 
offsets for each are given below. The hull radius at each station can be found by 
multiplying the offsets by half the maximum diameter, DU. 

If one were to use equations for true ellipsoids and parabolas, the entrance and the 
run would be too fine for a modem submarine. The displacement can be increased by 
using larger exponents («/and n a ), as in Equations (1) and (2). If even more displacement 
is required, a parallel middle body of cylindrical shape can be inserted at the maximum 
diameter. The prismatic coefficient, C p , is used to calculate volumes and for a cylinder the 
prismatic coefficient is 1. For a submarine-like body the prismatic coefficient can be 
evaluated in terms of its geometry. Using the above concept, the length of the parallel 
middle body (PMB) is the length overall (LOA) less 6 D, that is, LOA-6D. 
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Here x, and x„ are the distances from the maximum diameter. With these concepts, a very 
simple method of calculating the volume of the entire hull can be developed. This is true 
for the ends separately and for the PMB. Let V f , V „ , and Wmb denote, respectively, volume 
of the entrance, the run, and the parallel middle body, and let Q and C„ be the prismatic 
coefficients. The resulting equations are 

v , =^-( c „ 2AD ) ® 

V.-y(C,MD) (4) 
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The above can be combined into the following 



V = 



nD 



3.6C P a + — — 6 + 2AC P j 



(5) 



( 6 ) 



Figure 1 illustrates this concept [Ref. 4]. In this study our model submarine’s LOA is 
109.75 m. (360 ft.), D is 9. 15 m. (30 ft.), and the exponents n„ and n f are 3.0. 
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B. REGULAR WAVE RESPONSE 



In this study our main concern is the effects of surface waves on a near-surface 
submarine vertical motions. The assumptions here are that the fluid is ideal and the wave 




Figure 1 . Submarine Geometry [Ref. 4] 

and body motions are sufficiently small to linearize. Sea waves and ship motions are 
based on potential flow theory. In the simplest case it may be assumed that the waves 
incident upon the body are plane progressive waves of small amplitude, with sinusoidal 
time dependence. The solution to the water velocity distribution associated with the wave 
can be simplified to deep and shallow water approximations. When the wave length is 
greater than 20 times the depth, then shallow water approximations apply and depth 
becomes the controlling factor. The horizontal component of the velocity is not a function 
of depth, but is constant from top to bottom. The vertical component of water motion 
decays linearly from its maximum at the surface to zero at the bottom. The pressure under 
a shallow-water wave also is not a function of depth, but is just the hydrostatic pressure 
due to the amount of water above. When the water depth is greater than one quarter the 
wavelength then the deep water approximations apply and the water depth becomes 
unimportant. The horizontal and vertical components of velocity are equal and the orbits 
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become circles which decrease exponentially as a function of depth. Their motions 
become negligible at a depth equal to one half the wavelength [Ref. 5 ], Since we are 
interested in the vertical motion we consider the motions of a body that is allowed to 
heave and pitch only. Such motions are usually decoupled, for typical ships, from the 
horizontal plane motions in sway or yaw. The final coupled form of the heave and pitch 
equations of a ship in regular wave is then 

(m + ^33)773 + £35773 + C33773 + A35775 + £357/5 + C35775 = F 3 c'“ ( 7 ) 

(j. 55 + ^55)7)5 + £55775 + C55775 + A53773 + £53773 + C53773 = Fse ,a * ( 8 ) 

where m is the ship’s mass and I ss the mass moment of inertia with respect to the y axis. 
The Ajt terms correspond to added mass. The £,* terms correspond to hydrodynamic 
damping. Terms involving the coefficients ; i.e., C 33 , C 55 , and C S5 are related to tons per 
cm immersion, change in displacement per cm, and moment to trim one cm, respectively. 
The right hand side represents the heave, F 3 , Froude-Krylov and diffraction excitation 
forces. F 3 and F s are taken to be the complex exciting force and moment amplitudes, 
containing both amplitude and phase information. The previous equations of motion are 
valid for a ship with zero forward speed. If the ship possesses a forward speed U , this can 
be assumed, within linearity, constant. The only change in such a case is in the frequency 
0) due to a Doppler shift effect. In linear theory, the harmonic responses of the vessel, 
will be proportional to the amplitude of the exciting forces and at the same 
frequency, which is now ca (wave frequency of encounter) instead of 0). Consequently 
ship motions will have the form 
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T]j(t) = i(£kT].e im ,j = 3,5 (9) 

f}.(t) = -(0 2 e fj J e iM , 

where rj } is the complex response amplitude, and j= 3 for heave and j= 5 for pitch. 

Substituting equations (9) in (7) and (8), the e' m terms cancel out and the resulting 
equations are 

[~co 2 e (m + A 33 ) + i(o e B n + C , 3 - 4 ]r7 3 + (-m*A 35 + ico e B 35 + C 35 )t ] 5 = F 3 , (10) 

|^— tU f (/ 55 + A 55 } + i(D e B 55 + C 55 ji] 5 + (—G) e A 53 +ico e B 53 +C 53 )t ]3 = F 5 , (11) 

In equations (10) and (11) the origin is at the center of gravity, which is assumed to lie on 
the waterline. In the more general case, the term mf] 3 is substituted by mrj 3 -mx G f] 5 , 
and the term I 55 i) 5 by I 55 f] 5 +m(z G f] i -x G i) 3 ) , where 77, is the surge motion amplitude. 
In other words, a coordinate coupling is introduced. 

The determination of the coefficients and exciting forces and moments amplitudes 
represents the major problem in ship motions calculations. The problem can be simplified 
by applying a strip theory approach, where the ship is divided into transverse strips, or 
segments. The added mass and damping for each strip are relatively easily calculated, 
using two dimensional potential theory or by suitable two dimensional experiments. The 
sectional values are appropriately combined to yield values for A jk , B ik , C Jk , and F r 

To solve equations (10) and (11) for the complex amplitudes, the equations are 
written in the form 
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( 12 ) 



Pih+Qrj 5 = F„ 
Rrj 3 +Srj 5 = F S , 



where 



P = -co 2 e (m + A 33 ) + i(O e B n + C 3;t 

0 = Al5 **" ^35, 

/? — — <0,, A 53 + ico e B 5i + C 53 
S = — 0 ) e (/ 55 + A 55 ) + ico e B 55 + C 55 

The solutions to the coupled equations (12) and (13) is then given by 

- f,s-f 5 q 
_ f 5 p-f,r 

r\e = — — . 

5 PS-QR 



( 13 ) 



( 14 ) 

( 15 ) 



The ratio ( T|j /A ) is a quantity of fundamental significance, where A is wave amplitude, 
and is defined by Z/co, t/,0). Physically, this is the complex amplitude of body motion in 
the y-th mode, in response to an incident wave of unit amplitude, frequency to, and 
direction 0. The body itself moves with forward speed U. This ratio is generally known as 
the transfer function, or the response amplitude operator, RAO. The RAO can be 
calculated once the added mass, damping, exciting, and hydrostatic forces are known. 

The absolute vertical displacement at a point x along the length of the hull, due to 
heave and pitch is given by 

ZvA=*h-xri5’ ( 16 ) 
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and since t] 3 and rj 5 are the complex amplitudes in heave and pitch motion, respectively, 
^ VA contains both magnitude and phase information. Particularly of interest in this study 
is the relative vertical motion between a point in the ship and the surface of the 
encountered wave. The relative motion in regular waves is found by subtracting the free 
surface motion from the vertical ship motion at the desired point, taking account of their 
phase relationship. The free surface motion is composed of the incident wave, the 
diffracted wave, the radiated wave, and the Kelvin wave due to the ship’s steady forward 
speed. The traditional assumption is that the principal component is the incident wave; 
i.e., the incident wave is not distorted by the presence of the ship. Then the amplitude of 
the relative vertical motion in general is given by 

£*=«*- A*", (17) 

where A is the wave amplitude and k is the wave number. Then the RAO which requires 
only the scalar or absolute magnitude is 







*h~xri 5 Jkx 


A 




€ 

A 



The significance of the relative motion response is that the moments of their spectrum 
provide probability measures related to anticipated deck wetness, bow slamming or 
particularly for our study sail broaching and periscope submergence. 

Figures 2 and 3 show RAO amplitudes and phases versus wave to ship length 
ratio of our model submarine’s heave, pitch, and relative vertical motion at 3 submarine 
diameter depth and 5 knots forward speed when A equals 5 feet. 
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Amplitude (nondimensional) 




Figure 2. Amplitude of RAO for heave/pitch and relative motion 




Figure 3. Phase Angle of RAO for heave/pitch and relative motion 
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C. IRREGULAR WAVE RESPONSE 



Wave patterns in an open sea are ever changing with time and space, in a manner that 
appears to defy analysis be it linear or second order Stokes. Ambient waves on the surface 
of the sea are dispersive as well as random. Random refers to the character of the wave 
height distribution. The continuous distribution of sinusoidal waves have continuously 
distributed amplitude and phase so that in summation the variation of wave height with 
time is not systematic in any respect, but random. The practically useful data extractable 
from a random wave record h{t) is its spectral density, S(( o). The random h(t) record is 
processed in such a way to produce a curve of S( CO) versus wave frequency, (0. The 
spectral density is obtained from a wave height record taken over a time period for which 
the sea conditions are assumed to be unchanging, in an average sense (stationary). This 
corresponds to a certain sea state. The function 5(co,0) is called the spectral energy density 
or simply the energy spectrum. More specifically, this is a directional energy spectrum; it 
can be integrated over all wave directions to give the frequency spectrum 



Usually in the fields of ocean engineering and naval architecture it is customary to 
assume that the waves are long crested which means the fluid motion is two dimensional 
and the wave crests are parallel. With such a simplification it is possible to use existing 
information for the frequency spectrum (19), which is based on a combination of theory 
and full scale observations. 




( 19 ) 
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For most purposes we are interested primarily in the larger waves. The most common 
parameter that takes this into account is the significant wave height, H„ 3 , defined as the 
average of the highest one third of all waves. This can be computed from 






( 20 ) 



In this equation, mo is the area under the spectrum S(co) integrated over the entire range of 
frequencies CD. An average frequency of the spectrum can be defined as the expected 
number of zero upcrossings per unit time, that is, the number of times the wave amplitude 
passes through zero with positive slope. The final result here is 



CD Z = 



xl/2 

m 2 ' 
v m oy 



( 21 ) 



The average period between zero upcrossings is 



r = — = 2?r 

CD, 




( 22 ) 



More meaningful frequency parameters can be obtained from the set of moments, 
which depend on spectrum shape 



m n = J (O n S(a))dco , rc=0,l,2,... (23) 

o 

In particular, the area, m 0 , is the variance or the total energy of the spectrum. Also m 2 is 
variance of velocity and m 4 is variance of acceleration. 

A good model for fully developed seas is the classical Pierson-Moskowitz 
spectrum. This spectral form depends upon a single parameter which is the significant 
wave height. It is intended to represent point spectrum of a fully-developed sea. Fetch 
and duration are assumed to sufficiently large so that the sea has reached steady state, in a 
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statistical sense. This spectral family should be recognized as an asymptotic form, 
reached after an extended period of steady wind, with no contamination from an 
underlying swell. Using the spectral family, along with the similarity theory of S. A. 
Kitaigorodskii, Pierson and Moskowitz (1964) [Ref. 6] arrived at the following analytical 
formulation for ideal sea spectra. 



S, + (co) = one-sided incident wave spectrum 
g = acceleration of gravity 
Hvi = significant wave height 
CD = wave frequency 

In Figure 4 we can observe typical Pierson-Moskowitz wave spectra for 5 m. significant 




( 24 ) 



where 



wave height. 



50 



T 



45 - 




wave frequency (rad) 

Figure 4. Typical Pierson-Moskowitz wave spectra 



Any conclusions drawn on the seakeeping behavior of a ship based on the critical 
examination of motion response in regular waves can, at best, assume only academic 
significance. The establishment of the seakeeping behavior of a ship has to be done in a 
realistic seaway. With the spectral description of sea waves given before, we can return to 
the subject of body motions and generalize the results of regular harmonic waves. If the 
sea waves are described by the random distribution , and if the response of the body to 
each component wave is defined by a response amplitude operator Z(co,0), the body 
response will be 

77,(0 = *Jj Z.(m,£))e'“JA(m,6>) . (25) 

The principal assumption here is that linear superposition applies, as it must in any event 
for the underlying development of the RAO and the spectrum. 
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Like the waves themselves, the response (25) is a random variable. The statistic of 
the body response are identical to the wave statistics, except that the wave energy 
spectrum S is multiplied by the square of the RAO (this is a property of linear systems). 
Thus, if the subscript R represents any body response, we have 

Sr (®) = |Z R (<y)| 2 S (co) , (26) 

where Z*(co) is the RAO of the response R, and Sfco) the spectrum of the seaway. Equation 
(26) can then be utilized to obtain the spectrum of the response R. Figure 5 displays the 
spectrum of response of the relative vertical motion at the top of our model submarine’s 
sail while submarine’s forward speed is 5 Knots and it is at 3 submarine diameter depth. 
Also seaway is modeled by Pierson-Moskowitz spectrum with 5 m. significant wave 
height and head seas. 

To a large extent, equation (26) provides the justification for studying regular 
wave responses. The transfer function Z*(oo) is valid not only in regular waves, where it 
has been derived, but also in a superposition of regular waves, and ultimately in a 
spectrum of random waves. Generally speaking, a vessel with favorable response 
characteristics in regular waves will be good in irregular waves, and vice versa. 
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wave frequency (rad) 

Figure 5. Spectrum of response for relative vertical motion 
The average period between zero upcrossings was determined by equation (22), 
and the number between zero upcrossings per unit time is 



N« = — p_ 

2 2n\m£ ’ 



( 27 ) 



where m * , m * are the moments of the particular response R, whose spectral density is 
given by equation (26). Equation (27) can be generalized for the case of the average 
number of upcrossings above a specified level a as in 



n* =— h 

2 '“ 2k V m, 



J ex P 



a 



2 \ 



2mlj 



( 28 ) 



Equation (28) can be utilized to determine such events deck wetness and bow slamming 
for a surface ship or periscope submergence and sail broaching for a near surface 
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submarine. If/represents height of the periscope over calm sea surface level, the number 
of periscope submergence events per hour is 



where m 0 , m 2 are the moments of the vertical relative motion spectrum at periscope. The 
same equation can be used to estimate the frequency of sail broaching, with / substituted 
by the distance between top of the sail and encountered wave surface. Of course, m 0 , m 2 
are now the moments of the relative motion spectrum at top of the sail. When we 
calculate the frequency of vehicle collision with the sea-bed for shallow water operations 
this time m 0 , m 2 are the moments of the absolute motion spectrum at either the bow or the 
stem of the submarine and/ is its clearance from the sea-bed. 




( 29 ) 
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in. DEEP WATER RESULTS 



A. OPERABILITY INDEX 

As we mentioned in the first chapter in this study we considered two criteria when 
the submarine conducts near surface operations in deep water. One is the number of 
periscope submergence events per hour, which we can calculate with Equation (29) This 
number (N P i) 'was selected as 300, which corresponds to five periscope submergence 
events per minute. This is an arbitrary number and different numbers could be picked for 
different operational considerations. Since the same number was used for all cases, the 
results presented here are representative for all possible choices. The other criterion is the 
number of sail broachings per hour, which is also calculated from Equation (29). This 
number (N p2 ) was selected as one. This places far greater emphasis on sail broaching than 
periscope submergence since sail broaching is closely related to visual detection. 

Now we have the tools to compute the two performance indices defined above in 
a given seaway. Suppose that the submarine conducts periscope depth operation in a 
seaway characterized by a significant wave height (Pierson-Moskowitz Wave Spectrum), 
so that the sea spectrum is defined. For all round the clock boat headings relative to the 
predominant wave direction for which the operations are to be conducted, a polar plot 
diagram similar to the one in Figure 6 is prepared. Significant wave heights are 
represented along the radial direction of the polar plot. The shaded area in the plot shows 
wave height and wave direction combinations where the selected tactical assessment 
criterion is exceeded. Letting the polar area of the disk in Figure 6 be A 0 and the subset of 
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Ao within which the boat can conduct the operation be A, a performance index 
characterizing the ability of the boat perform this operation in the specified submarine 
velocity and depth can be defined as 100(AM o ). Generally a submarine’s forward speed 
ranges from three to twelve knots in periscope depth operations. We used three, five, 
eight and eleven knots submarine forward speeds, U, in our calculations. Depths, h, 
beneath the surface were selected from 2.5 to 5 boat diameters measured from the keel 
up. In computing the above index we could easily take into account the probability of 
occurrence of a particular sea state and wave heading angle in the area of interest by 
introducing appropriate weight factors. In this study we assumed that all possible sea- 
states and wave heading angles are equally probable. In the following sections we discuss 
the results for both criteria/operability indices and also for the combined 
criterion/operability index, where both criteria are taken into consideration at the same 
time. 




Figure 6. Typical performance assessment of a submarine 
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B. RESULTS OF PERISCOPE SUBMERGENCE CRITERION 

We begin by presenting results with regards to periscope submergence criterion 
alone. Typical polar plots are shown in Figures 7 through 30. The operability index is 
presented for different speeds and operating depths in Figures 3 1 through 40. Based on 
these results we can draw the following general conclusions: 

1. Head seas appear to result in a larger number of expected criterion violations than 
following seas. This is true regardless of the actual number used in establishing the 
criterion. 

2. For a given sea direction it is possible to reduce criterion violations for higher sea 
states. This simply means that the motion point moves more in phase with the 
incoming waves. It should be emphasized, however, that at such high sea states the 
average wave height may exceed the exposed periscope length. Since the periscope 
moves in phase with the waves, the operator’s visual horizon may be very small. This 
situation has been reported in practice and although the criterion is not exceeded, 
operations are very difficult to conduct. Such a situation can only be analyzed with 
proper visual simulation studies. 

3. An optimum operating depth can be found which minimizes the expected number of 
periscope submergence events. This depth depends on the forward speed, but it 
appears to be a weak function of speed. 
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Figure 7. Sea state-polar plot, showing SOE for a submarine {7=3 Knots, h=2.5D 
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Figure 8. Sea state-polar plot, showing SOE for a submarine U - 5 Knots, h=2.5D 
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Figure 9. Sea state-polar plot, showing SOE for a submarine U= 8 Knots, h=2.5D 
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Figure 10. Sea state-polar plot, showing SOE for a submarine U= 1 1 Knots, h=2.5D 
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Figure 1 1. Sea state-polar plot, showing SOE for a submarine U= 3 Knots, h=3D 
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Figure 12. Sea state-polar plot, showing SOE for a submarine U= 5 Knots, h=3D 
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Figure 13. Sea state-polar plot, showing SOE for a submarine U= 8 Knots, h=3D 
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Figure 14. Sea state-polar plot, showing SOE for a submarine U= 11 Knots, h=3D 
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Figure 16. Sea state-polar plot, showing SOE for a submarine U= 5 Knots, h=3.5D 



28 



120 



60 




90 8 




Figure 18. Sea state-polar plot, showing SOE for a submarine U= 1 1 Knots, h-3.5D 
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Figure 20. Sea state-polar plot, showing SOE for a submarine U= 5 Knots, h=4D 
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Figure 22. Sea state-polar plot, showing SOE for a submarine U=\ 1 Knots, h=4D 
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Figure 23. Sea state-polar plot, showing SOE for a submarine U = 3 Knots, h=4.5D 
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Figure 24. Sea state-polar plot, showing SOE for a submarine U= 5 Knots, h=4.5D 
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Figure 25. Sea state-polar plot, showing SOE for a submarine U = 8 Knots, h=4.5D 
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Figure 26. Sea state-polar plot, showing SOE for a submarine U= 1 1 Knots, h=4.5D 
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Figure 27. Sea state-polar plot, showing SOE for a submarine U= 3 Knots, h=5D 
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Figure 29. Sea state-polar plot, showing SOE for a submarine U= 8 Knots, h=5D 
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Figure 30. Sea state-polar plot, showing SOE for a submarine U=l 1 Knots, h=5D 
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Figure 3 1 . 01 vs. Submarine Speed Plot at 2.5 D depth 




Figure 32. 01 vs. Submarine Speed Plot at 3 D depth 
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Figure 33. 01 vs. Submarine Speed Plot at 3.5 D depth 




Figure 34. 01 vs. Submarine Speed Plot at AD depth 
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Figure 35. 01 vs. Submarine Speed Plot at 4.5 D depth 




Figure 36. 01 vs. Submarine Speed Plot at 5 D depth 
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Figure 37. OI vs. Submarine Depth Plot in Submarine Diameters at 3 Knots 




Figure 38. OI vs. Submarine Depth Plot in Submarine Diameters at 5 Knots 
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Figure 39. 01 vs. Submarine Depth Plot in Submarine Diameters at 8 Knots 




Figure 40. 01 vs. Submarine Depth Plot in Submarine Diameters at 1 1 Knots 
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C. RESULTS OF SAIL BROACHING CRITERION 



Typical polar plots for the sail broaching criterion event are shown in Figures 41 
through 56. The operability index is presented for different speeds and operating depths 
in Figures 57 through 65. Based on these results we can draw the following general 
conclusions: 

1. There does not appear to be a consistent dependence of the operability index on sea 
direction. Certain sea directions, however, greatly reduce the operability index for 
some speed/depth combinations. 

2. For a given sea direction, higher sea states correspond to smaller operability indices. 

3. The operability index does not vary significantly with speed or depth, and it generally 
increases with increasing depth.. 
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Figure 41. Sea state-polar plot, showing SOE for a submarine U= 3 Knots, h=2.5D 
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Figure 42. Sea state-polar plot, showing SOE for a submarine U= 5 Knots, h=2.5D 
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Figure 45. Sea state-polar plot, showing SOE for a submarine U = 3 Knots, h=3D 
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Figure 46. Sea state-polar plot, showing SOE for a submarine U= 5 Knots, h=3D 
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Figure 47. Sea state-polar plot, showing SOE for a submarine U= 8 Knots, h=3D 
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Figure 48. Sea state-polar plot, showing SOE for a submarine f/=l 1 Knots, h=3D 
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Figure 49. Sea state-polar plot, showing SOE for a submarine U= 8 Knots, h=3.5D 
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Figure 50. Sea state-polar plot, showing SOE for a submarine U= 1 1 Knots, h=3.5D 
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Figure 52. Sea state-polar plot, showing SOE for a submarine U = 5 Knots, h=4D 
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Figure 53. Sea state-polar plot, showing SOE for a submarine U = 8 Knots, h=4D 
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Figure 56. Sea state-polar plot, showing SOE for a submarine (7=1 1 Knots, h=4.5D 




Figure 57. 01 vs. Submarine Speed Plot at 2.5 D depth 
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Figure 58. 01 vs. Submarine Speed Plot at 3D depth 




Figure 59. 01 vs. Submarine Speed Plot at 3.5D depth 
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Figure 60. 01 vs. Submarine Speed Plot at AD depth 




Figure 61. 01 vs. Submarine Speed Plot at 4.5D depth 
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62. 01 vs. Submarine Depth in Submarine Diameters at 3 Knots 




Figure 63. 01 vs. Submarine Depth in Submarine Diameters at 5 Knots 
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64. 01 vs. Submarine Depth in Submarine Diameters at 8 Knots 




Figure 65. 01 vs. Submarine Depth in Submarine Diameters at 1 1 Knots 
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D. RESULTS OF COMBINED CRITERIA 



Results for both criteria combined are presented in Figures 66 through 99. The 
individual conclusions that were drawn previously continue to hold in this case as well. It 
appears that certain selections for speed and operating depth may result in higher values 
for the operability index. It should be mentioned, however, that this depends on the 
relative magnitude of the individual criteria, N p i and N P 2 . As an example, consider Figure 
97. The overall operability index at 5 knots and depth 4.5 diameters appears to be less 
than 10%. As Figure 83 demonstrates, however, this is entirely due to the periscope 
submergence criterion violations. If this were of no concern in this case, the operability 
index would be raised to 100% as shown in Figure 61. 




270 

Figure 66. Sea state-polar plot, showing SOE for a submarine U - 3 Knots, h-2.5D 
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Figure 68. Sea state-polar plot, showing SOE for a submarine U= 8 Knots, h=2.5D 
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Figure 69. Sea state-polar plot, showing SOE for a submarine U= 11 Knots, h=2.5D 
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Figure 7 1 . Sea state-polar plot, showing SOE for a submarine (7=5 Knots, h=3D 
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Figure 73. Sea state-polar plot, showing SOE for a submarine U= 1 1 Knots, h=3D 
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Figure 74. Sea state-polar plot, showing SOE for a submarine U=3 Knots, h=3.5D 
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Figure 75. Sea state-polar plot, showing SOE for a submarine U= 5 Knots, h=3.5D 
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Figure 76. Sea state-polar plot, showing SOE for a submarine U= 8 Knots, h=3.5D 



59 



60 




270 



Figure 77. Sea state-polar plot, showing SOE for a submarine U= 11 Knots, h=3.5D 
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Figure 78. Sea state-polar plot, showing SOE for a submarine U - 3 Knots, h=4D 
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Figure 79. Sea state-polar plot, showing SOE for a submarine U=5 Knots, h=4D 
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Figure 80. Sea state-polar plot, showing SOE for a submarine C7=8 Knots, h=4D 
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Figure 81. Sea state-polar plot, showing SOE for a submarine U-\ 1 Knots, h=4D 
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Figure 83. Sea state-polar plot, showing SOE for a submarine (7=5 Knots, /i=4.5£> 
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Figure 84. Sea state-polar plot, showing SOE for a submarine U = 8 Knots, h=4.5D 
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Figure 85. Sea state-polar plot, showing SOE for a submarine U= 1 1 Knots, h=4.5D 
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Figure 86. Sea state-polar plot, showing SOE for a submarine (7=3 Knots, h=5D 
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Figure 87. Sea state-polar plot, showing SOE for a submarine U= 5 Knots, h=5D 
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Figure 88. Sea state-polar plot, showing SOE for a submarine <7=8 Knots, h=5D 
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Figure 89. Sea state-polar plot, showing SOE for a submarine U= 1 1 Knots, h=5D 




Figure 90. 01 vs. Submarine Speed Plot at 2.5 D depth 
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Figure 91. 01 vs. Submarine Speed Plot at 3 D depth 




Figure 92. 01 vs. Submarine Speed Plot at 3.5D depth 
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Figure 93. 01 vs. Submarine Speed Plot at 4 D depth 




Figure 94. 01 vs. Submarine Speed Plot at 4.5£> depth 
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Figure 95. 01 vs. Submarine Speed Plot at 5 D depth 




Figure 96. OI vs. Submarine Depth in Submarine Diameters at 3 Knots 
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97. 01 vs. Submarine Depth in Submarine Diameters at 5 Knots 




Figure 98. OI vs. Submarine Depth in Submarine Diameters at 8 Knots 
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Figure 99 . 01 vs. Submarine Depth in Submarine Diameters at 1 1 Knots 
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IV. SHALLOW WATER RESULTS 



In this chapter we present results for shallow water operations; i.e., for finite water 
depth beneath the keel. Shallow water corrections are implemented in the software 
package as described in [Ref. 3]. The results are presented in terms of the two criteria 
used in deep water, periscope submergence and sail broaching as well as for an additional 
criterion which quantifies the probability of collision between the boat and the sea bed. 
This is calculated based on the absolute vertical motion at either the bow or the stem and 
the clearance from the sea bed at the appropriate motion point. All results are presented 
for a standard depth between keel and free surface of three diameters. 

A. RESULTS OF PERISCOPE SUBMERGENCE CRITERION 

Figures 100 through 108 present the results for the periscope submergence 
criterion. Each polar plot is for a given depth beneath the keel and for four forward 
speeds, namely three, five, eight, and eleven knots. Figure 109 shows the operability 
index in percentage normalized with respect to its deep water value, as a function of both 
speed and shallow depth. The percent change in the operability index is calculated by 
subtracting the shallow water value from the deep water value and dividing by the deep 
water index. Therefore, positive values show that the operability in shallow water is 
smaller compared to deep water. It can be seen that the operability index does not change 
much and it does not follow a consistent trend. For certain speeds it decreases in shallow 
water, while for other speeds, it increases and then starts to decrease. The numerical 
values are quite small and the shape of the polar plots does not change significantly for 
different water depths. It appears, therefore, that the operability index with regards to 
periscope submergence remains essentially the same for deep and shallow waters. 
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Figure 100. Sea state-polar plot, showing SOE for different U’s and 5 D Water Depth 




Figure 101. Sea state-polar plot, showing SOE for different (7‘s and 6 D Water Depth 
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Figure 102. Sea state-polar plot, showing SOE for different U's and ID Water Depth 




Figure 103. Sea state-polar plot, showing SOE for different U’s and 8 D Water Depth 
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Figure 104. Sea state-polar plot, showing SOE for different U’s and 10Z) Water Depth 




Figure 105. Sea state-polar plot, showing SOE for different U’s and 15D Water Depth 
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Figure 106. Sea state-polar plot, showing SOE for different U ' s and 20 D Water Depth 




Figure 107. Sea state-polar plot, showing SOE for different U ' s and 30 D Water Depth 
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Figure 108. Sea state-polar plot, showing SOE for different U ' s and 40D Water Depth 




Figure 109. Change in 01 with respect to infinite depth vs. Water Depth 
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B. RESULTS OF SAIL BROACHING CRITERION 



Results for the sail broaching criterion are presented in Figures 1 10 through 118. 
The same parameters as for periscope submergence were used here as well. The percent 
change in the operability index is shown in Figure 1 19. Here it seems that water depth has 
a more important effect on both the value of the operability index and the shape of the 
polar plots. In general, the operability index is decreasing as the water depth becomes 
smaller. 




Figure 110. Sea state-polar plot, showing SOE for different t/’s and 5 D Water Depth 
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Figure 111. Sea state-polar plot, showing SOE for different Vs and 6 D Water Depth 




Figure 1 12. Sea state-polar plot, showing SOE for different U's and ID Water Depth 
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Figure 113. Sea state-polar plot, showing SOE for different U’s and 8 D Water Depth 




Figure 1 14. Sea state-polar plot, showing SOE for different U’s and 10Z) Water Depth 
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Figure 115. Sea state-polar plot, showing SOE for different U ’ s and 15D Water Depth 




Figure 1 16. Sea state-polar plot, showing SOE for different (/’ s and 20D Water Depth 
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Figure 1 17. Sea state-polar plot, showing SOE for different C/’s and 30 D Water Depth 




Figure 1 18. Sea state-polar plot, showing SOE for different £/’ s and 40D Water Depth 
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Figure 1 19. Change in 01 with respect to infinite depth vs. Water Depth 



C. RESULTS OF COLLISION CRITERION 

Results for the sea-bed collision criterion are summarized in Figures 120 through 
125. The operability index is shown in Figure 126. As expected, the operability index 
decreases with decreasing water depth. The shape of the polar plots, i.e., specific sea state 
and sea direction combinations, also varies significantly with water depth. 
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Figure 120. Sea state-polar plot, showing SOE for different U's and 5 D Water Depth 




Figure 121. Sea state-polar plot, showing SOE for different U ' s and 6D Water Depth 
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Figure 122. Sea state-polar plot, showing SOE for different f/’s and ID Water Depth 




Figure 123. Sea state-polar plot, showing SOE for different U’s and 8 D Water Depth 
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Figure 124. Sea state-polar plot, showing SOE for different U’s and 10D Water Depth 




Figure 125. Sea state-polar plot, showing SOE for different If s and 15D Water Depth 
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Figure 126. 01 vs. Water Depth 



D. RESULTS OF COMBINED CRITERIA 

Results for all three criteria combined are presented in Figures 127 through 162. 
Lighter shades correspond to the periscope submergence criterion while the darkest 
shades correspond to the collision criterion. It can be seen that, in general, the sail 
broaching criterion dominates the collision criterion for the parameters selected; i.e., one 
occurrence per hour for sail broaching and one occurrence per two hours for collision. 
Should different values for the parameters had been selected, the collision criterion might 
dominate instead. The total change in the operability index is shown in Figure 163. It can 
be seen that, in general, the operability index is decreasing for decreasing water depth. 
Shallow water effects seem to be insignificant for depths exceeding 30 submarine 
diameters. 
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Figure 128. Sea state-polar plot, showing SOE for U - 5 Knots and 5 D Water Depth 
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Figure 129. Sea state-polar plot, showing SOE for U = 8 Knots and 5 D Water Depth 
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Figure 130. Sea state-polar plot, showing SOE for U= 11 Knots and 5 D Water Depth 
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Figure 131. Sea state-polar plot, showing SOE for (7=3 Knots and 6D Water Depth 
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Figure 132. Sea state-polar plot, showing SOE for U= 5 Knots and 6 D Water Depth 
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Figure 134. Sea state-polar plot, showing SOE for (7=1 1 Knots and 6 D Water Depth 



92 




60 



90 8 




270 

Figure 136. Sea state-polar plot, showing SOE for U= 5 Knots and ID Water Depth 
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Figure 137. Sea state-polar plot, showing SOE for U= 8 Knots and ID Water Depth 
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Figure 138. Sea state-polar plot, showing SOE for U= 1 1 Knots and ID Water Depth 
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Figure 139. Sea state-polar plot, showing SOE for t/=3 Knots and 8 D Water Depth 
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Figure 140. Sea state-polar plot, showing SOE for U=5 Knots and 8 D Water Depth 
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Figure 141. Sea state-polar plot, showing SOE for (7=8 Knots and 8D Water Depth 
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Figure 142. Sea state-polar plot, showing SOE for U= 11 Knots and 8 D Water Depth 
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Figure 143. Sea state-polar plot, showing SOE for U= 3 Knots and 10Z) Water Depth 
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Figure 144. Sea state-polar plot, showing SOE for £7=5 Knots and 10£> Water Depth 
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Figure 145. Sea state-polar plot, showing SOE for U= 8 Knots and 10 D Water Depth 
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Figure 146. Sea state-polar plot, showing SOE for U=l 1 Knots and 10D Water Depth 
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Figure 147. Sea state-polar plot, showing SOE for U= 3 Knots and 15D Water Depth 
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Figure 148. Sea state-polar plot, showing SOE for U = 5 Knots and 15D Water Depth 
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Figure 149. Sea state-polar plot, showing SOE for (7=8 Knots and 15D Water Depth 
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Figure 150. Sea state-polar plot, showing SOE for (7=1 1 Knots and 15D Water Depth 
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Figure 151. Sea state-polar plot, showing SOE for f/=3 Knots and 20 D Water Depth 
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Figure 152. Sea state-polar plot, showing SOE for U = 5 Knots and 20D Water Depth 
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Figure 154. Sea state-polar plot, showing SOE for (7=1 1 Knots and 20 D Water Depth 
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Figure 155. Sea state-polar plot, showing SOE for (7=3 Knots and 30 D Water Depth 
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Figure 156. Sea state-polar plot, showing SOE for U= 5 Knots and 30D Water Depth 
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Figure 157. Sea state-polar plot, showing SOE for 17=8 Knots and 30 D Water Depth 
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Figure 158. Sea state-polar plot, showing SOE for U= 11 Knots and 30D Water Depth 
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Figure 159. Sea state-polar plot, showing SOE for U= 3 Knots and 40 D Water Depth 
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Figure 160. Sea state-polar plot, showing SOE for U=5 Knots and 40D Water Depth 



105 



120 



60 



90 8 




270 



Figure 162. Sea state-polar plot, showing SOE for U= 1 1 Knots and 40D Water Depth 
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Figure 163. Change in 01 with respect to infinite depth vs. Water Depth 
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V. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

In this study we have evaluated the vertical plane response of submersible 
vehicles in the proximity of a free surface in deep and shallow waters. A potential flow, 
strip theory solver is used for the computations. For deep water operations periscope 
submergence and sail broaching criteria are used to map the operability envelope of the 
vehicle for different sea directions, sea states, operating depths and vehicle speeds. The 
operability of the vehicle is quantified by calculating a certain operability index. For 
shallow water operations in addition to the above two criteria, a third criterion is 
considered which is the collision of the vehicle with the sea-bed. The primary conclusions 
from this study are summarized below: 

1 . For the periscope submergence criterion head seas appear to cause larger number of 
violations than following seas regardless of the water depth. Also, the operability 
index does not change much with water depth. The effect of shallow water on this 
criterion is insignificant. An optimum operating depth can be found which minimizes 
the expected number of periscope submergence events. This depth is a weak function 
of vehicle speed. 

2. For the sail broaching criterion the operability index does not appear to depend on sea 
direction in a consistent way. Higher sea states correspond to smaller operability 
indices for a given sea direction. The operability index does not change significantly 
with speed or operating depth and it generally increases with increasing operating 
depth. For the sail broaching criterion, water depth has a more important effect than 
for periscope submergence on both the value of the operability index and the shape of 
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the polar plots In general, the operability index is decreasing with decreasing water 
depth. 

3. For the collision criterion in shallow water, the operability index decreases as the 
water depth becomes smaller. The shape of the polar plots changes significantly with 
water depth. 

4. For all criteria combined, it appears that certain combinations of vehicle speed and 
operating depth may result in higher values for the operability index. It should be 
mentioned that this depends on the relative magnitude of the individual criteria. In 
general, the sail broaching criterion dominates the collision criterion for the 
parameters selected in this study. It can be seen that, in general, the operability index 
is decreasing for decreasing water depth. Shallow water effects seem to be 
insignificant for depths exceeding 30 submarine diameters. 

B. RECOMMENDATIONS 

The following is a list of recommendations for further research on near surface 

response of submersible vehicles:' 

1. For periscope submergence in a given sea direction the motion point appears to move 
more in phase with the incoming waves as sea states become more severe. Even 
though the criterion is not exceeded in such high sea states, the average wave height 
may exceed the exposed periscope length. Since the periscope moves in phase with 
the waves, the operator’s visual horizon may be very small. This can cause the 
operations to be difficult to conduct even though the criterion is not violated. Such 
situations should be analyzed with proper simulation studies. 
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2. Evaluate the effect of different geometric hull parameters, on the operability index. 
Such parameters may be, for example, the length, diameter, and prismatic coefficients 
of the hull. 

3. Evaluate the effects of second order wave forces and motions on vehicle response. 
Even though these motions are slowly varying and can be controlled to a certain 
extent, they may alter both the values of the operability indices and the shape of the 
corresponding polar plots. 



ill 
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